Fermionic collective excitations in a lattice gas of Rydberg atoms 
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We investigate the many-body quantum states of a laser-driven gas of Rydberg atoms confined to 
a large spacing ring lattice. If the laser driving is much stronger than the van-der-Waals interaction 
among the Rydberg sates, these many-body states are collective fermionic excitations. The first 
excited state is a spin-wave that extends over the entire lattice. We demonstrate that our system 
permits to study fermions in the presence of disorder although no external atomic motion takes 
place. We analyze how this disorder influences the excitation properties of the fermionic states. Our 
work shows a route towards the creation of complex many-particle states with atoms in lattices. 
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Experiments with ultracold ground state atoms have 
opened a doorway to the study of many-particle systems 
and revealed deep insights into many-body physics and 
the dynamics of phase transitions [l| . Entangled states of 
atomic ensembles serve as a resource for quantum infor- 
mation processing, precision measurement and the gener- 
ation and storage of light 2] . Recently, there is a growing 
interest in highly excited atomic (Rydberg) states as the 
strong state-dependent interaction between these atoms 
allows the implementation of fast quantum information 
protocols [1, HI, [j| and the creation of strongly interacting 
(spin) systems which show collective excitation behavior 
[fj,l3]- The long-ranged character of the interaction allows 
to entangle atoms with a spatial separation of the order of 
several micrometers, as has recently been demonstrated 
experimentally d,H|. 

Motivated by this development we env isag e atoms 
trapped in the sites of a deep large spacing [h], H3 ring 
lattice with lattice constant a ~ [im (see Fig. QJi). We 
consider bosonic atoms occupying the ground state of 
the respective potential well. Rydberg states are excited 
by a laser. The typical timescale associated with the fast 
(electronic) Rydberg dynamics is of the order of hundreds 
of nanoseconds. In this regime the external motion of 
the atoms is frozen and the time-evolution of a Rydberg 
gas is described by an arrangement of spin 1/2 particles 
with Ising type interaction [lj, [TH, 0, EH . Here the two 
spin states correspond to the atomic ground and Rydberg 
state, and the laser detuning and its Rabi frequency give 
rise to an effective magnetic field. 

Here we show that the above-mentioned setup allows 
to create and explore collective many-particle states that 
are entangled and extend over the entire lattice, e.g. spin 
waves, on a microsecond timescale. We consider a sce- 
nario in which the interaction with the laser is stronger 
than the interatomic interaction. This strong driving 
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FIG. 1: a: Ring lattice with spacing a being much larger than 
the extension a of the Wannier functions (deep lattice). The 
internal atomic degrees of freedom at each site are described 
by the (collective) states \P) k and \R) k - b: Level structure 
in the regime in which Q*, = 2> P and |A| <C fi. The 
spectrum splits into manifolds which can be labeled by the 
quantum number m of the operator £\ ai . For sufficiently 
large O, the coupling between manifolds established by Hi 
and H2 can be neglected. The (constrained) dynamics inside 
the m-subspaces is then determined by -ff xy - 



constrains the system's dynamics such that the collec- 
tive excitations of the atomic ensemble are described by 
spinlcss fermions. We characterize these many-particle 
states and discuss how they can be excited from an ini- 
tial state in which only ground state atoms are present. 
We finally show that the system offers the possibility to 
study fermions in the presence of a disorder potential al- 
though no external atomic motion takes place. Here the 
disorder is created by a spatially randomly varying Rabi 
frequency. We focus on a situation where this variation is 
caused by fluctuations of the atom number in each lattice 
site. These fluctuations are generated by quenching a su- 
pcrfluid gas of ground state atoms by a sudden increase 
of the lattice depth. 

The Hamiltonian of a frozen Rydberg gas in a ring 
lattice with L sites is given by 
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with Vk = (1/2) (1 + (ri^) being the projector onto the 

(k) 

state in which a Rydberg atom is present on site k. a x 

and cri are the Pauli spin matrices and ai L+1 ^ — cri. 
The first two terms of Eq. |T]) account for the cou- 
pling of the atoms to the excitation laser within the 
rotating-wave approximation. If more than one atom is 
present per lattice site a simultaneous excitation of two 
or more Rydberg atoms is forbidden due to the large 
(interaction induced) energy shift of this multiply ex- 
cited state. This implies that in general, on each site 
k, the excitation laser couples the product state \P) k = 
[Ifiv/Ji ® ■■■ ® OflWjv to tne superatom state \R) k = 

^"* 1/2 «5{[l»->*]i ® [l3> fe ] 2 ® - ® l\9) k ] Nh } (S is the sym- 
metrization operator, N k is the number of atoms con- 
tained in site k and \g) and \r) are the sing le atom ground 
and Rydberg states, respectively) [S Eff. Throughout, 
we consider these two states to constitute a (super)atom 
located at site k, and we define a { z k) \P) k = - \P) k and 

CT 2 fc ' l-^)fc = \R)k- They are coupled by the collective 
Rabi frequency S1& = y/Nkflo, where Slo is the single- 
atom Rabi frequency. The detuning of the laser with re- 
spect to the transition \P) k — > \R) k is A. The last term 
of Hamiltonian (JTJ) describes the interaction between (su- 
per)atoms located at different sites. We consider a regime 
in which the lattice spacing a is much larger than the ex- 
tension of the Wannier functions a (see Fig. QJ,). The in- 
teraction strength between (super) atoms in neighboring 
sites is then given by (3 = Ce /a 6 , with Cq being the van- 
der-Waals coefficient. We consider only nearest neighbor 
interaction, which is well justified for sufficiently large 
rings since the next nearest neighbor interaction is a fac- 
tor of 64cos 6 (7r/L) smaller. Hamiltonians of the form 
(P) have been extensively and successf ully used to s tudy 
strongly interacting Rydberg gases [ISESESESESEif- 
Throughout this work we are interested in many- 
particle states that emerge in the limit of strong laser 
driving, i.e Q k 3> (3 and |A| <C Sl k . Since in this 
regime the first term of Eq. fl} dominates, it is con- 
venient to diagonalize it applying the unitary transfor- 
mation U = Y[ k exp(— *( 7r /4)cy fe ' 1 ). We then find H = 
WH spin U = H xy +H 1 +H 2 + (3L/A with 
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where H* y is the famous xy-model of a spin chain 
with a transverse magnetic field [III]. After the uni- 
tary transformation, the eigenstates of ai are - in 
terms of the (super)atom states - given by |±) fc = 

(l/V2)W [\P) k ± \R) k ] with a[ k) \±) k = ± |±) fc . 



Let us first consider the limit where each site of the 
lattice is occupied by the same number of atoms, i.e., 
N k = No, which implies a constant Rabi frequency 
fife = SI. This is achieved, for example, if the lattice 
is initialized in a Mott-insulator state. For SI ^> (3, the 
spectrum of H decays into manifolds of nearly degener- 
ate states having the same eigenvalue m with respect to 
the operator ai k K These manifolds are separated by 
approximately 20 (see Fig. [T|d). As indicated in Fig. 
[TJj, coupling between the manifolds is only caused by Hi 
and H 2 . The corresponding transition rates are approx- 
imately A 2 /SI and (3 2 /fl for Hi and H 2 , respectively. 
Hence, for sufficiently strong driving the system's dy- 
namics is constrained to the m-manifolds and Hi 2 can 
be neglected. 

The dynamics is then governed by the icy-model Hamil- 
tonian ([2|) which can be solved by introducing creation 
and annihilation operators of spinless fermions (c k and 
Cfe, respectively) through the Jordan- Wigner transforma- 
tion c fc = exp (i^Y^jZi [111- One obtains 
J? xy = Hq + Hi, with 
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and H h = -f (c\c x - c L c[\ (e l7m + + l) being the 
boundary term due to the ring. The latter depends 
on the operator n + — X)fe=i c I c fc ano - counts the num- 
ber of fermions / (super)atoms in the state |+). The 
creation operators rjj, n of the eigenexcitations of Hamil- 
tonian ([l} can be found by a Fourier transform. Due 
to the boundary term one has to distinguish between 
an odd (p = 6) and even (p = e) number of fermions: 

vhn = 7T Efe=i ex P ( ia n k ) 4 with a n = 2nn/L and 
a c n — 2{n — l/2)ir/L. The diagonal Hamiltonian reads 
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Its ground state reads \G) = Y\ k \ —} k and the excited 
states are of the form \ij) = vliVlj\G) or \ijk) — 

vl.ivljvl k 1^)- T nese states are eigenstates of the ini- 
tial Hamiltonian fl} in the limit Q k = £1 3> (3. 

We now address the two questions: ' Which states can 
be excited from a given initial state?'' and How can they 
be excited?'. To this end we define the initial conditions: 
i) The system shall be in a state in which no Rydberg 
atoms are present, i.e., |0) = Y[ k \P)k- ^) The laser is 
turned off, i.e., Slo(O) = and set to a certain detuning 
A(0) = Ao. This corresponds to the left-hand side of the 
plots in Fig. ® where we present the energy spectrum of 
Hamiltonian |T]) (L = 10) as a function of SI. 

Which states can be excited? - The symmetry proper- 
ties of the system impose certain selection rules which im- 
ply that not all fermionic excitations are accessible from 
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the initial state |0) by a sweep of the laser parameters. 
In Ref. [15| it was shown that Hamiltonian (Q]) and also 

(|6|) are invariant under cyclic shifts (X^a^X — crj fc+1 ^) 

and reversal of the lattice sites {Ttfa^lZ = cr^ L k+1 '). 
Since the initial state is also invariant under the action 
of both operators, i.e., X< |0) = TV |0) = |0), only states 
from this (fully symmetric) subspace are accessible when 
A and f2o are varied. 

Let us first study which single fermion states are ac- 
tually accessible. A general state is given by r]l n \G) 
where n — 1 . . . L. For a fully symmetric state we require 
vl,JG) = OmjG) = Oh^OO^G) = 0^l n O\G) 
where O is a placeholder for X and 1Z. Here we have 
used the fact that \G) - just like |0) - is a fully sym- 
metric state. Direct calculation shows that Ttfrfc n 1Z = 

™+ and X^r]l n X = e-^rjl^^ 1 + 
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-j=c\(e l7rn+ — 1). The only fully symmetric state contain- 



ing one fermion is thus |1) = F° r states that 

contain two fermions one finds |2 n ) = r/l L _ n+1 r)\ n \G) 
where n = 1 . . . [L /2J . In order to give a physical mean- 
ing to the states it is more instructive to use operators 
that refer to the atomic basis: 
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Hence, |1) is a spin wave, or a superatom state that 
extends over the entire lattice and the states |2„) are 
entangled two-atom states. Such excitations are of in- 
terest as they can yield a resource for the creation 
of single photons and photon pairs [20|. Excitations 
containing three fermions are constructed via \3i pq ) — 

75 l r )li r lip r )l q ~ vIl-ivIl^pVIl^I where the in- 
dices I, p and q can assume all distinct combinations 
obeying I + p + q = eL with e = 1,2. Further excited 
states are constructed in a similar fashion. A thorough 
investigation of their properties will be presented else- 
where. 

How can the many-particle states be excited? - We seek 
to selectively populate states in the regime fl ^> (5. To 
this end we vary the laser parameters as a function of 
time, i.e., A — > A(t) and n — > Slo(t). Our starting point 
are the above-mentioned initial conditions (left-hand side 
of Fig. [2]). We now consider two possible excitation meth- 
ods which depend on the initial sign of the detuning: 

In case of Ao < (left inset of Fig. [2]) the ground 
state at f2 = does not coincide with the initial state |0). 
With increasing fl the energy of this state cuts through a 
number of avoided crossings. It is thus not a trivial task 
to find a proper trajectory (A(£), fio(t)) which eventually 
connects |0) with a desired eigenstate of H xy . Only in 
certain cases (see Refs. [13, [3) such a sequence can be 
guessed. A more general framework for finding a proper 




FIG. 2: Energy spectrum (fully symmetric states) of Hamilto- 
nian for a lattice with 10 sites in units of (3. Many-particle 
states can be accessed from the state |0) by a temporal vari- 
ation of the laser parameters. If Ao > (background figure 
and right inset [magnified view]) the state |0) is adiabatically 
connected to the ground state |G) of H xy (thick orange line). 
See text for further explanation. If Ao < (left inset) the 
state |0) cuts through a number of avoided crossings as Q 
is increased. The dashed line shows a passage to the state 

| 523456 }• 



trajectory is provided by Optimal Control theory [2l| . 
Here certain constraints on the trajectory can be imposed 
(e.g. to avoid too fast temporal oscillations) and the 
desired fidelity with which the final state is achieved can 
be set. This method is successfully applied to quantum 
information processing [2^ |. molecular state preparation 
[23| and optimization of number squee zing of an atomic 
gas confined to a double well potential [24j |. 

We pursue a different route and set Ao > at fi = 0. 
As shown in the right inset of Fig. fJl |0) now coin- 
cides with the ground state. Ramping up f2 this state 
is adiabatically connected to the ground state \G) of 
H xy . For small the distance to adjacent states can 
be increased by increasing A such that non-adiabatic 
transitions are strongly suppressed. From the state \G) 
many-particle states can be excited by introducing an 
oscillating detuning A(t) = A osc cos(ojt) giving rise to a 
time-dependent H±. Within the rotating- wave approx- 
imation the coupling matrix element between \G) and 
the fully symmetric single fermion state is (G\ Hi |1) = 
— (A osc V / L)/4. The oscillating detuning acts like a 'laser' 
which drives transitions between different many-particle 
states. Hence, if the (radiofrequency) oscillation of the 
detuning is on resonance, i.e., u> « ujl = 2£! + (3/2 (see 
Fig. [2]), Rabi flops between the ground state and the 
first excited state are performed. Resonance with the 
|1) — » |2„) transitions can be avoided by choosing a not 
too high value of O. Here, level shifts due to H2, which 
are of the order /3 2 /fi, bring such unwanted transitions 
out of resonance. This offers also a possibility to excite 
higher lying states by first performing a 7r-pulse on the 
\G) — > |1) transition followed by a 7r-pulse with frequency 
adapted to the transition |1) — ► \2 n ). 

Due to the limited lifetime of the Rydberg atoms (e.g. 
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FIG. 3: Absorption profile (excitation probability of a single 
fermion state) for the \G) — > |1) transition for two strengths 
of disorder, Q//3 = 10 and L = 50. The disorder strength is 
controlled by the mean number of atoms No per site. The 
results are averages over 1000 realizations. 



66 us for Rb in the 60s state), the whole excitation proce- 
dure has to take place on a (is time scale. Once prepared, 
the many-particle states are mapped to a stable state by 
performing a global 7r-pulse on the single-atom transition 
|?') — > |s), where |s) is a stable storage state. Experimen- 
tally, the excited many-particle states can be detected 
by counting the number of atoms in the |+) -state which 
corresponds to the number of fermions. 

So far we have assumed a constant Q k . We will now 
consider a situation in which £l k is not constant but fluc- 
tuates randomly around a mean value, i.e. fl k = £l + Sfl k . 
In Eq. ([5]) the fluctuating part 6£l k introduces a ran- 
dom single particle potential for the fermions, and give s 
rise to the Hamiltonian of Anderson localization [251 ] . 
Hence, a lattice gas of Rydberg atoms offers the pos- 
sibility to study fermions in a disorder potential al- 
though no external atomic motion takes place. The 
required spatial variation of Q k can, for instance, be 
achieved by a speckle potential or standing waves with 
incommensurate frequencies [26l [27| . We present an al- 
ternative route: We start in a situation in which N s 
ground state atoms are prepared in a superfiuid state 

|SF) = {NglL^/ 2 )- 1 Efc=i 6 fe l vac ) where b l cre- 
ates a ground state atom in state \g) at site fc. Disor- 



der is introduced by a quench of this superfiuid through 
a sudden increase of the depth of the lattice potential. 
For a sufficiently large number of atoms N g the Rabi 
frequency of site k is then given by Q k = fioV-^fc — 
v^^o + Sn k = O + Sfl k with Sfl k = SN k n/(2N ) 
and No = N s /L. Here SQ k is a random variable with 
(Sfl k ) = and (SQ k 5Q m ) = fl 2 /(4N ) [S km - L~ x ] . The 
disorder destroys the symmetry properties of the system 
and hence also the selection rules for transitions between 
many-particle states. As a consequence, there are now in 
general L single fermion states (instead of 1) accessible 
from the state \G) when an oscillating detuning is ap- 
plied. Instead of a single sharp line at lot, the (averaged) 
absorption profile of the \G) — ► |l)-transition broadens 
and becomes asymmetric, see Fig. [3J For small disorder, 
one finds from second order perturbation theory, that the 
red wing is approximated by I(u>) ~ il 2 / (4No\uj — lol\ 2 )- 
This disorder-induced line broadening can be detected 
by counting the number of |+) -atoms as a function of 
the excitation frequency uj. 

In our considerations we have assumed that the atoms 
are strongly localized, i.e., o > <r. In practice there is 
a finite width of the wave-packet, caused by the uncer- 
tainty principle and finite temperature. This will lead to 
disorder also in the interaction energy /3, which can be 
also treated in the present framework but is beyond the 
scope of this work. Moreover, the number of atoms in a 
single site cannot be made too large due to three-body 
losses. Bearing this in mind, the system opens exciting 
perspectives for creating complex many-particle states 
with interesting prospects for the study of disorder and 
the generation of non-classical light. 
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